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Abstract 

The catalytic activity of cytochrome c (cyt c) to peroxidize cardiolipin to its oxidized form is required for the release of pro- 
apoptotic factors from mitochondria, and for execution of the subsequent apoptotic steps. However, the structural basis for 
this peroxidation reaction remains unclear. In this paper, we determined the three-dimensional NIVIR solution structure of 
yeast cyt c Y67H variant with high peroxidase activity, which is almost similar to that of its native form. The structure reveals 
that the hydrogen bond between MetSO and residue 67 is disrupted. This change destabilizes the sixth coordination bond 
between heme Fe^* ion and MetSO sulfur atom in the Y67H variant, and further makes it more easily be broken at low pH 
conditions. The steady-state studies indicate that the Y67H variant has the highest peroxidase activities when pH condition 
is between 4.0 and 5.2. Finally, a mechanism is suggested for the peroxidation of cardiolipin catalyzed by the Y67H variant, 
where the residue His67 acts as a distal histidine, its protonation facilitates 0-0 bond cleavage of H2O2 by functioning as an 
acidic catalyst. 
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Introduction 

Recently, it has been proposed that proteins can inherently 
possess a variety of conformations in vivo, and hence have 
fianctional diversity [1]. However, it's difficult to understand the 
mechanisms involved in diiTerent conformers from native protein 
to "non-native" protein [2-3], because it's hard to obtain the 
single conformer of the protein at different states. Cytochrome c 
{i.e., cyt c), a well-known electron transfer hemoprotein, has 
several conformers to carry out its biological functions beyond 
respiration [4] . Two of them have given rise to special attentions. 
One is the so-called "alkaline conformer" with weak peroxidase 
activity [5-6], in which the axial ligand of the heme iron, Met80, is 
replaced by Lys72, Lys73 or Lys79 at alkaline pH condition. The 
structural investigation on the alkaline conformer implied that this 
conformer might work as an electron transfer gate, and a folding 
intermediate [7] . The other is the "pro-apoptotic" conformer with 
enhanced peroxidase activity [8-10]. It can catalyze peroxidation 
of cardiolipin (a mitochondria specific phospholid) to its oxidized 
form, which is essential for the release of pro-apoptotic factors 
from mitochondria, and for performing the subsequent apoptotic 
processes. But the structural basis of this catalytic peroxidation 
remains unclear. 

In cyt c, Tyr67 is a highly conserved residue except that in the 
mitochondrial cyt c-558 of Euglena gracilis, where the residue 67 



is phenylalanine [11]. The residue Tyr67 is extremely important in 
the electron transfer because its side-chain extends into the heme 
pocket and forms a hydrogen bond network with axial ligand 
Met80, conserved residues Asn52 and Thr78 and a buried water 
molecule identified in the X-ray structure of cyt c (pdb code 
2YCC) [12-13], as shown in Figure 1. This hydrogen bond 
network not only helps to maintain the sulfur atom of Met80 in its 
heme iron ligand position and thus controls the redox potential of 
the protein, but also modulates the flexibility of the nearby peptide 
chain segment [i.e., residues 65 to 72), therefore affects chemical 
properties of the protein. More interestingly, Tyr67 is located at 
the distal position of heme pocket, the distance between Tyr67 - 
OH group and heme Fe '^ is 4.3 A, which is close to the distances 
of 5.84 A and 5.55 A between N£2 atom of the distal histidine and 
heme iron ion in horseradish peroxidase (i.e., HRP, pdb code 
1H5A) and cytochrome c peroxidase [i.e., CcP, pdb code 2CYP) 
[14-15], respectively. Thus, Tyr67 was suggested to be a possible 
apoptotic trigger via a conformation resulted by site-directed 
mutagenesis [16-17]. 

Originally, in order to convert cyt c into a peroxidase, we 
successfully designed cyt c Y67H and Y67R variants [16-17]. 
Replacement of Tyr67 by His67 is to introduce a distal histidine, 
while mutation from Tyr67 to Arg67 is to produce a distal 
arginine. The results from electronic spectrum (a weak charge 
transfer band centered at 695 nm for Fe-S bond) and CD 
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Figure 1. Structure of the heme region of cyt cand the location 
of Tyr67 in the structure (pdb code: 1YIC), showing the 
hydrogen bond network formed by amino acids His18, 
Tyr67, MetSO, Asn52 and Thr78. 

doi:10.1371/journal.pone.0107305.g001 

spectrum (loss of a negative band centered at 4 1 6 nm) implied that 
the heme Fe-S (MetSO) bond was not broken in the Y67H variant, 
but really ruptured in the Y67R variant. Moreover, the midpoint 
temperatures (T^J of the Y67H variant (50°C) and of the Y67R 
variant (30°C) are much lower than those of the wild-type cyt c 
(higher than 80°C) and of its Y67F variant (90°C) [18], indicating 
that the detachment of Met80 to heme iron is significantly easier in 
the Y67H variant than that in wild-type cyt c. Both Y67H and 
Y67R variants have high peroxidase activities [16-17]. Thus, the 
Y67H and Y67R mutants stand for a typical pro-apoptotic 
conformer of cyt c. Since there is no solid structural evidences to 
interpret the basis of peroxidation in this state, it's necessary to 
study their three dimensional conformation. Thus, in this paper, 
by conventional two-dimensional H- H NMR methods [19], we 
determined the NMR solution structure of the cyt c Y67H variant. 

In this paper we investigate the structure of the Y67H mutant. 
This mutant was constructed using the pBTR2 plasmid which 
encodes the yeast iso-l-cyt c (CYCl) genes. This is a cytochrome c 
variant carrying two point mutations K72A/C102T. The change 
from Cysl02 to Thrl02 is used to prevent disulfide dimerization of 
the protein and protein auto-reduction [7,20-2 1] , which may alter 
the cytochrome's thermodynamic and spectroscopic parameters. 
The mutation from Lys72 to Ala72 is utihzed to prevent this 
residue from serving as a ligand in the alkaline form of cyt c at 
pH 7.0. For the sake of simplicity, throughout this paper, the 
double-site K72A/C102T variant is used as the reference and 
referred to as native yeast iso-1 cyt c. The triple mutants Y67H/ 
K72A/C102T and Y67R/K72A/C102T used in tiiis study are 
referred to the Y67H and Y67R variants, respectively. 

Results 

Heme destruction upon native cyt c and its Y67H variant 
reacting with H2O2 

To better understand the mechanism of peroxidation by the 
Y67H variant, we first probed the stability of native cyt c and its 
Y67H variant only in the presence of H2O2. As shown in Figure 2 
(For clear description, the concentration of H2O2 in Figure 2B is 



1 mM, because the absorption at 410 nm decays extremely fast 
when it is 10 mM), in the absence of guaiacol, once H2O2 reacts 
with native cyt c and its Y67H mutant, the visible absorption 
maximal values at ~408 nm (Soret band), 530 nm (a+fi bands) 
and ~353 nm (5 band) become weaker and weaker, indicating 
that the heme ring in both proteins is destructed very quickly [22]. 
This observation is apparently different from myoglobin [22-23], 
which in the absence of substrate for peroxidation reacts with 
H2O2 to generate compound II as confirmed by the absorption 
band at ~420 nm (compound II structure is shown in Figure 3). 
Figure 4A demonstrates that the native cyt c and its Y67H mutant 
have different behaviors vs pH variation. For the Y67H variant, 
the observed rate (A„i,s) of heme destruction is the highest at about 
pH 4.0. When pH value is equal to or higher than 6.0, the ko\K 
values for the Y67H variant are ahnost kept at 0.02 s~', but is 
dramatically decreased when pH is less than 4.0. For the native cyt 
c, the A„|,s value of heme destruction is not completely changed 
upon pH condition being varied. Our previous acid titration 
studies on both native cyt c and its Y67H variant [1 7] implied that, 
diflFerent from wild-type cyt c, at pH 4.0-5.2 (especially at pH 4.0), 
Fe-S bond was clearly cleaved since the absorbance at 695 nm 
underwent a sharp decrease. The cleavage of Fe-S bond generates 
a vacancy for H2O2 coordination, so that heme ring is easier to be 
oxidized by H2O2. Thus, the A„bs value of heme destruction of the 
Y67H variant at pH between 4.0 and 5.2 is significantiy higher 
than those at other pH conditions. 




300 350 400 450 500 550 600 650 
Wavelength (nm) 



Figure 2. Heme destruction of the wild-type cyt c (A) and its 
Y67H variant (B) as indicated by dissipation of the Soret band. 

Spectra were scanned every 20 sec, and the arrows indicate the 
direction of change. The concentrations of H2O2 were 1 0 mM in (A) and 
1 mM in (B), the pH condition in both cases is 6.0. 
doi:1 0.1 371/journal.pone.01 07305.g002 
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Figure 3. The possible mechanism for guaiacoi reaction with H2O2 catalyzed by the cyt c Y67H mutant at pH 4.0. 

doi:l 0.1 371/journal.pone.01 07305.g003 



Steady-State kinetic studies 

Aft(;r we performed heme destruction by H2O2, we then 
investigated the peroxidase activity of wild-type cyt c and its Y67H 
mutant on guaiacoi vs pH change. As shown in Figure 4B, upon 
pH being varied, the native cyt c displays no significant changes in 
the peroxidase activity, but the Y67H mutant demonstrates the 
high(;st peroxidase activity when pH is between 4.0 and 5.2, 
increased by about 4—10 folds, compared to that of native cyt c. As 
we know, at pH 4.0-5.2 the axial Fe-S bond is cleaved, generating 
a channel for loading H2O2 and substrate guaiacoi which 
enhances the peroxidation reaction. 

Sequence-specific assignment 

Using two dimensional (2D) 'H-'H total correlation spectros- 
copy (TOCSY) and correlation spectroscopy (COSY) in H2O and 
D2O for spin patterns and 2D 'H-'H nuclear Overhauser effect 

spectro.scopy (NOES"Y) spectra in H2O for sequential NH-NH and 
Ha-HN connectivities, we made the assignments for the Y67H 
mutant under our experimental conditions (20°C and 50 mM 
phosphate buffer, pH 7.0), according to the extensive lists of 
assignments for yeast cyt c reported in the literatures [24—26]. The 
chemical shifts of residues Glnl6, Thrl9, Gly29, AsnSl and Leu32 
were firsdy identified because their spin patterns were resolved 
outside the diamagnetic envelope. The sequence-specific assign- 
ment was successfully performed in the regions of residues 3-13, 
50-54, 61-68, 71-74 and 88-102, referencing with the data for 
the native cyt c [24—26]. Residues, especially close to the heme 
ring and the mutation site, exhibit different chemical shifts, which 
were finally assigned and given in Table 1. In total, more than 
88% of the expected proton resonances had been assigned for the 
Y67H mutant. 



NMR signals assignment of axial ligands and heme of the 
Y67H variant 

One dimensional (ID) 'H NMR spectra shown in Figure 5 
indicate that the hyperfine shifts of the Y67H mutant are 
completely distinct from those of the native cyt c, so NMR signals 
of heme of the Y67H variant cannot be assigned directly according 
to the similarity of the ID 'H-NMR spectra between the native cyt 
c and its Y67H mutant. The complete assignments of hyperfine 
shifted signals of the heme protons, the axial ligands His 18 and 
Met80, and the heme bound Cysl4 and Cysl7 of Y67H were 
achieved through analysis of the 2D 'H-'H NOESY spectra 
tailored to the relaxation properties of these paramagnetic 
resonances and through one-dimensional nuclear Overhauser 
effect (NOE) experiments, as reported before [24—25]. The 
characteristic peaks at —31.0 ppm (peak 3b) and —22.7 ppm 
(peak 2b) in the high-field shifted region of ID ' H-NMR spectrum 
of native cyt c (Figure 5) were previously assigned as the axial 
ligand MetSO Hyl and s-CHj [24-25], respectively. Their 
corresponding peaks are at —27.8 ppm (peak 3a) and — 
19.5 ppm (peak 2a) in the high-field shifted region of ID 'H- 
NMR spectrum of die Y67H mutant. Thus, NMR data suggest 
that the Met80 is still served as the sixth axial ligand in Y67H 
mutant. 

The well resolved signals at 33.1, 14.0, 28.4, and 17.7 ppm 
(Figure 5B) were suggested to be the four methyl groups of the 
heme of Y67H mutant. The signal at 14.0 ppm can be assigned to 
I-CH3 because it has NOEs with the signals of a thioether group 
and methyl groups of Leu32, Leu68 and Leu94, which were also 
<)bs(;r\ed in native cyt c. The resonances at 33.1 and 14.0 ppm 
display strong NOE connectivities to a signal at —0.3 ppm. The 
only proton which is equidistant from two methyl groups is the 
meso-H5, positioned between the I-CH3 and the 8-CH3 groups. 
The assignment of the signal at 33.1 ppm as 8-CH3 is consistent 
with its connectivities with 52-CH3 of Leu32 and S-CH3 of Met64. 
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Figure 4. Steady-state kinetics studies on (A) lieme destruction 
and (B) peroxidase activity vs phH change of wild-type cyt cand 
its Y67H mutant. 

doi:1 0.1 371/journal.pone.01 07305.g004 

The remaining resonances are assigned as 3-CH3 and 5-CH3. The 
resonance at 17.7 ppm was assigned as 5-CH3 due to its NOE 
connectivities to Gly29 Ha (at -0.5 ppm) and Pro30 H5 (at — 
1.6 ppm). The resonances of the other heme substituents were 
assigned by analyzing the NOESY connectivities, starting from 



position 1 and 3 and "hopping" around the macrocycle (Table 2 
and Figure SI in File SI). 

The peaks at the positions of —26.32 ppm (peak lb) and 
24.54 ppm (peak 6b) (Figure 5) which belong to the side-chain 
aromatic protons He 1 and H(52 of the axial ligand His 1 8 in native 
cyt c [24—25], respectively, are shifted to the positions (peak la, — 
23.85 ppm; peak 6a, 22.65 ppm) in ID 'H-NMR spectra of the 
Y67H mutant (Figure 5). The saturation of the signal at — 
23.85 ppm allows the detection of a NOE with a proton 
resonating at 11.10 ppm, which displays NOE connectivities to 
two H;S protons of His 18 resonating at 15.55 ppm and 8.36 ppm. 
The assignments of two H;8 protons of His 18 are consistent with 
the NOE patterns between them and the amide HN protons of 
Leu32 and Thrl9 observed in native cyt c [24—25]. So the signals 
at 11.81 ppm and —23.85 ppm were assigned as the side chain 
atoms H(5l and Hsl of Hisl8, respectively. The detectabUity of the 
signal H(51 of His 18 indicates that its H-bond with Pro30 is 
preserved in the Y67H mutant. The NOE patterns between the 
two H/? protons of His 18 and the signal at 22.65 ppm imply that 
this signal is from H(52 of His 18, which also displays a NOE 
connectivity with a thioether group at 2-position (1.47 ppm) in the 
heme ring. The saturation of the signal at 11.10 ppm indicates 
that it has NOE patterns with two H/J protons of Hisl8, suggesting 
that it comes from the backbone NH of His 1 8, which is further 
confirmed by the NOE peaks between it and the Ha of Cysl7 and 
Thrl9 HN. The assignment of Ha of Hisl8 at 8.84 ppm was 
performed based on NOE patterns between this signal and the HN 
proton, H/i!2 proton of His 18 and between this signal and HN of 
Thrl9. The ID 'H-NMR spectra of the Y67H mutant show 
broad signals at —19.5 ppm (peak 2a in Figure 5A) and 
22.65 ppm (peak 6a in Figure 5B), both do not exchange with 
deuterium in D2O solution. The line widths of these signals suggest 
that they are from protons of the metal ligands, most probably the 
axial ligands. Thus, based on the assignments of native cyt c [24— 
26], these were appointed as Met80 E-CH3 and His 18 H52, 
respectively. 

NMR solution structure determination of the Y67H 
mutant 

The helical conformations of the Y67H mutant were charac- 
terized by using the reported techniques [24—26], further ensured 
by using observed NOE patterns (Figure S2 in File SI), are 
displayed in the regions of 3-13, 50-54, 61-68, 71-74, and 88- 





1 


02 aa. They are almost similar to those presented in the solution 


Table 1. 


New assignments of Y67H mutant. 






residues 


chemical shift (ppm) 




His'" 


HN 11.10, Ha 8.84, H^l 8.36, Hj92 15.55, Hi51 11.81, Ht)2 22.65, Hi:l 


-23.85 


Pro^° 


Ha 4.12, Hyl -0.71, H-,'2 0.01, H(5 -1.61 




Leu" 


<51-CH3 0.62, (52-CH3 1.08 




lie" 


y-CHj -0.12, (5-CH3 0.23, y-CHj 0.49, 1.03 




Trp=' 


He 7.52 






HN 8.30, Ha 428, Hj31 3.12, Hj32 3.18, H(5 6.82, HC2 8.00 






Ha 3.36, H/J1 0.23, H/J2 1.05, H-;' 0.68, 15I-CH3 -3.05, 152-CH3 -1.02 




Met™ 


HN 8.67, Hyl -27.8, E-CH3 -19.5 




Phe'^^ 


HN 8.44, Ha 4.56, H/jl 2.98, H/j2 3.25 




Leu^^ 


HN 8.16, Ha 4.25, H/il 0.72, H/i2 1.13, Wy 0.83, lil-CHj 0.19, 152-CH3 


0.54 



doi:l 0.1 371 /journal.pone.01 07305.t001 
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Figure 5. ID ^H-NWIR spectra of cyt c. (A) The high-field shifted region of the native cyt c (down) and its Y67H variant (upper) spectra. (B) The 
down-field shifted region of the native cyt c (down) and its Y67H variant (upper) spectra. The peaks in ID ^H-NIVIR spectrum were assigned as (la) 
His18 HeI, (2a) His71 Hel, (3a) His71 H62, (4a) heme 8-CH3, (5a) His18 H82, (6a) heme 3-CH3, (7a) heme 5-CH3, (8a) His18 H|32, (9a) heme I-CH3 and 
(1 Oa) heme 7-Ha2, respectively. The peaks in 1 D ^ H-NIVIR spectrum of native cyt c were assigned as (1 b) Hisi 8 HeI , (2b) MetSO E-CH3, (3b) Met80 Hy, 
(4b) heme 8-CH3, (5b) heme 3-CH3, (6b) His18 H52, (7b) heme l-Hal, (8b) His18 HP2, (9b) heme 7-Ha1 and (10b) His18 H81, respectively. 
doi:1 0.1 371/journal.pone.01 07305.g005 

structures of the wild-type Saccharomyces cerevisiae iso-1 cyt c [24— 
26] . The helix IV (7 1-74 aa) in native cyt c is disappeared in some 
conformers of the Y67H mutant. 

The solution structures of the Y67H variant were determined by 
using meaningful and acceptable 1699 NOEs through the 
program XPLOR [27]. The heme group, the axial ligands, and 
the two cysteine (Cysl4 and Cysl7) covalently attached to the 
porphyrin are treated as new patch residues. Among them, the 
patch residue HEC stands for cyt c-type heme ring, PHEM is 
served as one axial ligand histidiiie of heme iron, and PHMT 
represents methionine worked as another axial ligand of heme 
iron. PHCB and PHCC are used as the residues Cysl4 and Cysl7. 
These applications for solution structure calculation of cyt c-type 
heme contained proteins are integrated in the XPLOR-NIH 
software package. One hundred structures were initially calculat- 
ed. Upon no NOE violations more than 0.3 A being found, 30 
pairs of hydrogen bond constraints involving slowly exchanging 
amide protons consistently presented in the initially calculated 
structures, were introduced as further constraints in the final stage 
of the structure calculations. Finally, the best 20 structures with the 
lowest energy selected from 100 calculated structures, were used to 
represent the three-dimensional structures of the Y67H mutant. 
The conformers of this bundle show no NOE violation more than 
0.3 A, and have backbone atoms root mean square deviation 
(RMSD) values of 0.77±0.15 A, heavy atoms RMSD values of 
1.15 ±0.17 A with respect to their corresponding mean structures 
(calculated for residues —5-103 and heme group), respectively. 
Figure 6 demonstrated the backbone superimposition of 20 
structures and ribbon representation of three-dimensional solution 
structures of the Y67H mutant. The breakdown of experimental 
constraints per residue was summarized in Figure S3 in File SI. 
The solution structure infers the secondary structures predicted 
from NOE pattern mentioned above (Figure S2 in File SI). The 
entire structure statistics for the 20 conformers were summarized 



Table 2. The assignment of heme protons of cyt c Y67H 
mutant. 




atoms 


cyt c y67H (ppm) 


I-CH3 


14.0 


2-CH 


1.47 


2-CH3 


-1.25 


meso- H„ 


3.96 


3-CH3 


28.4 


4-CH 


3.02 


4-CH3 


2.38 


meso- Hji 


-0.39 


5-CH3 


17.7 


6- Hal 


-1.05 


6- Ha.2 


2.42 


6- H/il 


1.05 


6- H(32 


2.35 


meso- H,. 


7.71 


7- Hal 


11.09 


7- Ha2 


13.18 


7- H/il 


0.51 


7- H/i2 


1.30 


8-CH3 


33.1 


meso- H,i 


-0.31 



doi:1 0.1 371 /iournal.pone.Ol 07305.t002 
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in Table 3, where 84.4% of the residues were located in tlie most- 
favored regions of the Ramachandran plot for the Y67H mutant. 
These data indicate that the solution structure of the Y67H 
mutant is reasonable. 

Discussion 

Mutation from Tyr67 to His67 changes hydrogen bond 
network in heme pocket 

To investigate the conformational changes in heme moiety 
relative to heme ring produced by the mutation from Tyr67 to 
His67, in Figure 7, the superrmpositions were performed by 
overlaying backbone Co; atoms in secondary structural regions and 
heme ring including all carbon and nitrogen atoms to form the 7i- 
conjugate porphyrin system. The conformers of native cyt c (pdb 
code lYIC) and its Y67H mutant are almost identical with a 
global backbone atoms RMSD value of 1 . 1 5 A, which can explain 
the fact the exposure of heme group of the Y67H variant to the 
solvent is quite similar to that of the native form (only 0. 1 % 
increment in the Y67H variant evaluated with the program 
MOLMOL, when considering the iron ion and the carbon and 
nitrogen atoms to form the 7i-conjugate porphyrin system.). In 
native cyt c, -OH group of Tyr67 forms one hydrogen bond to 
Met80 sulfur, holding an internal molecule of water. In the Y67H 
mutant, although the side-chain of His67 is almost in the same 
position of Tyr67 side-chain (Figure 8), but the distance between 
His67 Ns2 and Met80 sulfur is close to 4.0 A in aU fmal 20 best 
conformers (obviously larger than that between oxygen atom of 
Tyr67 -OH and sulfur atom of Met80 is 3.34 A in the native cyt c) 
(Table SI). Moreover, the distance between His67 N£2 and heme 
coordinated Fe ion was measured as 5.48 A, much longer than 
that (4.30 A) between the oxygen atom of Tyr67 -OH group and 
Fe'^ ion in the native cyt c; the distances in the Y67H mutant 
between Thr78 Oyl atom and oxygen atoms Olot and 02a of 
heme 6-propionate are 5. 1 1 A and 6.07 A, respectively, both 
larger than those (4.55 A and 3.44 A) in the native cyt c. The 
distances between Asn52 N52 and oxygen atoms Ola and 02a of 
heme 7-propionate were measured as 5.51 A and 3.87 A in the 
Y67H mutant, respectively, also larger than those (2.86 A and 
4.38 A) in the native cyt c. AU these differences in distances 
indicate that the mutation from Tyr67 to His67 disrupts the 



hydrogen bond between Met80 and residue 67, and also changes 
the hydrogen bond network. These structural changes weaken the 
axial Fe-S bond (the Fe-S bond length is 2. 1 9 A in native cyt c, and 
2.43 A in the Y67H mutant), the stability of the Y67H mutant is 
thus decreased. The weakened Fe-S bond results in slight shift of 
helix IV (71-74 aa) away from the heme ring (Figure 7A). This 
accords with the fact that the distance between Lys73 atom 
and Fe'^"^ ion is 19.6 A in the Y67H mutant, close to that (19.0 A) 
in the Y67F mutant [18], but clearly longer than that (16.9 A) in 
native cyt c. Therefore, our NMR structure of the Y67H variant 
reveals a basis why Y67H is more unstable than its wild-type. 

To probe the structural difference between the Y67H and Y67F 
mutants, we superimposed their conformations (Figure 7B) and 
measured the distance changes between some key residues 
(Figure 8B). The Y67H variant has a global backbone atoms 
RMSD value of 1.49 A with the Y67F variant (pdb code: ICTY), 
implying that their global foldings are almost identical. In 
comparison, helix IV (71-74 aa) only slightly shifts away from 
the heme pocket in the Y67H variant, which makes the heme 
cavity a litde looser than that in the Y67F variant. This is 
confirmed by the following observations. The mutation from 
Tyr67 to Phe67 destroys the hydrogen bond between Met80 and 
heme Fe''"^ (the bond length of Fe-S is 2.44 A, close to that in the 
Y67H variant), but the protein stability is increased by forming a 
more compact conformation than that in native cyt c [7]. 
Moreover, the distances between Trp59 Nsl and heme 7- 
propionate carboxyl group Ola atom was measured as 4.46 A 
in the Y67H variant, larger than that (3.86 A, hydrogen bond is 
observed) in the Y67F variant. The distance between Thr78 Oyl 
atom and the carboxyl atom 025 in heme 6-propionate is 6.07 A 
in the Y67H variant, much larger than that (3.08 A) in the Y67F 
mutant. These changes imply that the hydrogen bonds observed in 
the Y67F variant between Thr78 and heme 6-propionate, and 
between Trp59 and heme 7-propionate may disappear in the 
Y67H variant. 

As we knew, from the native cyt c {i.e. the conformer III, pdb 
code lYCC or lYIC) to its alkaline form (also called as conformer 
IV, pdb code ILMS) [7,24-26], the conformation of the sixth 
ligand is switched from Met80 in the native form to Lys73 in the 
alkaline form, a pH dependent conformational toggling, the 
ligation between side-chain NHj of Lys73 and heme iron results in 




Figure 6. 3D solution structures of cyt £rY67H variant. (A) Backbone of the family of 20 Y67H structures shown as a tube of variable radius. The 
radius is proportional to the RMSD of each residue. (B) Ribbon diagram representation of the Y67H mutant. The figures are generated with the 
program MOLIVIOL [48]. The iron ion was represented as ball in red. 
doi:1 0.1 371/journal.pone.01 07305.g006 
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Table 3. Experimental restraints and structural statistics for cyt c Y67H variant. 




parameters 


20 structures of the Y67H variant 


Distance restraints from NOEs 


Total NOE 


1699 


Intra-residue (i-j = 0) 


642 


Sequential (|i-j| = 1) 


392 


Medium range (1<|i-j|<5) 


213 


Long range (|i-j|>5} 


452 


H-bond pairs restraints 


30 


Structural statistics 


r.m.s.d versus tiie mean structure(A) 


All backbone atoms 


0.77±0.15 


All heavy atoms 


1.15±0.17 


Backbone atoms (secondary structure} 


0.48 ±0.1 6 


Heavy atoms (secondary structure) 


0.93±0.17 


r.m.s.d from the experimental restraints 


NOE distances (A) 


0.0193±0.00099 


RMSD from idealized geometry 


Bonds(A) 


0.00167±0.000055 


Angles(') 


0.23±0.0a72 


lmpropers(°) 


0.36±0.0088 


Ramachandran analysis 


Residues in most favored regions 


84.4% 


Residues in additionally allowed regions 


11.1% 


Residues in generously allowed regions 


3.3% 


Residues in disallowed regions 


1.1 


Number of bad contacts/100 residue^ 


0 


Overall G-factoi^ 


0.18 


"The programs PROCHECK and PROCHECK-NMR were used to check the overal 
''For the PROCHECK statistic, less than 10 bad contacts per 100 residues, and 
doi:1 0.1 371 /journal.pone.01 07305.t003 


quality of the structure and GLY and Pro are excluded from the Ramachandran analysis, 
an overall G-factor larger than —0.5 are expected for a good quality structure. 




Figure 7. Upon overlaying Ca atoms in second structural region and heme backbone atoms, the conformational comparison: (A) 
between wild-type cyt c(grey, pdb code 1YIC] and its Y67H variant (green); (B) between the Y67F (magenta, pdb code 1CTY) and 
Y67H (green) mutants; (C) between Y67H (green) and alkaline form (red, pdb code 1LMS). 

doi:1 0.1 371 /journal.pone.01 07305.g007 
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B 



K73^ 





Figure 8. The position comparison of some l<ey residues (including Asn52, Trp59, Tyr67, Lys73, Tlir78, Phe82) relative to lieme ring 
upon overlaying heme ring: (A) between the native cyt c(grey, pdb code 1YiC) and its Y67H variant {green); (B) between the cyt c 
Y67H variant {green) and the Y67F variant {magenta, pdb code 1CTY); (C) between cyt cY67H variant {green) and the alkaline form 
of cyt c {red, pdb code 1 LMS). 

doi:1 0.1 371 /journal.pone.01 07305.g008 



the loss of helix IV (71-74 aa). In the Y67H mutant, Met80 is still 
ligated with Fe'^''", thus the alkaline form of cyt c and the Y67H 
mutant demonstrate big differences in the conformation with a 
global backbone atoms RMSD value of 2.64 A (calculated for 
residues 9-102 aa), as shown in Figxire 7C. The movement of the 
distal loop in the alkaline form is obviously larger than that in the 
Y67H mutant, the relative position of helix III (residues 63-69 aa) 
to the heme ring in Y67H mutant is closer than that in the alkaline 
form. The relative position between helix II (50-54 aa) and helix 
III in the Y67H mutant is also closer than that in the alkaline 
form, implying that the channel from heme pocket to the solvent 
in the alkaline form is open larger than that in the Y67H variant. 
This can be confirmed by the distance changes shown in 
Figure 8G and Table SI. The distances between Tyr48 -OH 
group and heme 6-propionate Olot and 02q( were measured as 
5.38 A and 7.59 A in the alkaline form, respectively, clearly larger 
than those (3.46 A and 5.29 A) in the Y67H variant. Similarly, the 
distances between Thr49 -OH group and heme 6-propionate 
Ola and O2o{ were measured as 7.26 A and 8.96 A in the alkaline 
form, respectively, much larger than those (3.68 A and 5.47 A) in 
the Y67H variant. Other residues, such as Trp59, Tyr67/His67, 
Thr78 also shift farther away from heme ring in the alkaline form 
than those in the Y67H variant. These structural difiFerences once 
again reveal that cyt c possesses functional diversity through 
conformational variability. 

Taken together, the current structure of the Y67H variant 
interprets why Y67H is the most unstable among native cyt c, its 
Y67H and Y67F mutants, and also the peroxidase activity of the 
Y67H variant, which wUl be discussed below. 

Mechanism of cardiolipin peroxidation catalyzed by the 
Y67H variant 

To successfully switch the classical electron transfer protein, cyt 
c, into a peroxidase, firstly we need to understand the main 
differences on the structural elements between cyt c and 



peroxidase, which were listed below: (1) the 6' axial ligand of 
heme in peroxidase is vacant for substrate loading; (2) in the heme 
pocket of peroxidase there is a distal histidine facilitating the 
formation of compound I intermediate (Figure 3); (3) in the heme 
pocket of peroxidase the existence of a distal arginine to assist this 
process. Tyr67 is located at the distal position of heme pocket of 
the yeast iso-1 cyt c, thus, we thought that the mutations from 
Tyr67 to His67, or to Arg67, might induce a distal histidine or a 
distal arginine, which could enhance peroxidation activity of cyt c. 
Within our anticipation, in guaiacol peroxidation at pH 7.0 and 
16°C, compared to the wild type cyt c, the catalytic efficiency 
i<:„, for the Y67H variant (where A,,„/^m = 89.38 M"'s"', 
i<:„, = 0.0096 M, A„„ = 0.858 s"') is increased by about 16 folds 
(for wild-type cyt c, the value k^^f/K^-^ is 5.12 M~'s~', 
= 0.0127 M, = 0.065 s"'), while W-^m of the Y67R 
variant is raised by about 8 times (for the Y67R variant, the value 



37.47 M" 



= 0.0162 M, A„, = 0.607 s"') [16- 



17]. The KJK^ value for the Y67H variant at 42°C (A,^„,/ 



K„, = 984.83 M" 



= 0.0087 M, L 



3.568 s" 



about 10-fold greater than that at 16°C, while for the Y67R 
variant, the k^^/K^ value (A„t/^m= 1844.44 M~'s~', 
i<:„, = 0.0090 M, k,^^= 16.600 s"') was almost 50-fold higher. In 
the Y67R variant, the mutation from Tyr67 to Arg67 may 
produce not only a distal arginine (although its structure is not 
available), but also a room for substrate coordination (because Fe- 
S bond is broken in the Y67R variant). Our current NMR 
structure of the Y67H variant (determined at pH 7.0) suggests that 
the sixth hgatioii bond Fe-S is attenuated due to the disruption of 
hydrogen bond between Met80 and Tyr67 existed in the native cyt 
c. The distance between Ns atom of His67 and heme ferric iron is 
measured as 5.48 A in the Y67H variant, very close to the 
distances of 5.84 A and 5.55 A between Ns2 atom of the distal 
histidine and heme iron ion in horseradish peroxidase (HRP, pdb 
codes 1H5A) and cytochrome c peroxidase (CcP, pdb codes 
2CYP)[14,15], respectively. Therefore, in the Y67H variant. 
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His67 may act as a distal histidine to increase the peroxidase 
activity. These results correlate with the peroxidase activity and 
conformational changes of cyt c resulted by the site-directed 
mutations under the physiological conditions (neutral pH) and 
without denaturant for the first time. 

However, the Y67H \ariant has much higher peroxidase 
activity at pH values from 4.0 to 5.2 than those at any other pH 
conditions (Figure 4), in contrast, native cyt c has very low 
peroxidase activity at any pH condition. How do we explain this 
phenomenon? In native cyt c, the hydrogen bond network, 
including the hydrogen bond between Tyr67 and Met80, stabilizes 
the ligation between MetSO and heme Fe'^'^ ion. Thus, the pH 
titration neither changes the hydrogen bond network, nor weakens 
the bond between MetSO sulfur atom and heme Fe''^ ion (this 
analysis is consistent with the results from UV-visible and CD 
spectroscopy [16-17]). While in the Y67H variant, the removal of 
the hydrogen bond between Tyr67 and Met80 weakens the heme- 
MetSO ligation, leading to that MetSO easily detaches from heme 
ring. The cleavage of Fe-S bond in the Y67H variant results in a 
channel in the heme crevice to H2O2 and substrate for 
peroxidation, which is a characteristic of the cardiolipin induced 
conformational transition in cyt c. Therefore, at pH values from 
4.0 to 5.2, two foUowing elements may stimulate the Y67H variant 
with the highest peroxidase activity: 1) His67 acts as a distal 
histidine; 2) vacancy generation specific for the substrate or H2O2 
loading. When pH is less than 4.0, the peroxidase activity of the 
Y67H variant is dramatically decreased due to the production of 
the low-spin intermediate His 1 S-heme-water, confirmed by a 
sharp increase at the absorbance at 620 nm, comparable to the 
wild- type cyt c [17]. 

At pH conditions from 4.0 to 5.2 in the absence of H2O2, the 
neutral aromatic ring of residue His67 of the Y67H variant is 
protonated (Figure 3). Thus, the protonated His67 in the Y67H 
variant direcdy acts an acid to cleave O-O bond in HgOj, losing a 
molecule of water. In contrast, the distal histidine in the 
peroxidases HRP or CcP is proposed to work as a dual acid- 
base catalyst in the peroxidation reaction (Figure S4 in File SI) 
[28], where the distal histidine at first functions as a general base to 
deprotonate H2O2, then protonated histidine acts as an acid to 
cleave O-O bond. By omitting the step to deprotonate H2O2, the 
protonation of His67 in the Y67H variant facilitates the 
peroxidation reaction of H2O2 with heme Fe ^ ion. Figure 3 
displays a mechanism about how the compound I is formed 
through the reaction catalyzed by the Y67H variant at pH 4.0. In 
the absence of the substrate guaiacol for peroxidation, the 
compound I reacts with H2O2, resulting in quick heme ring 
destruction. Otherwise, the compound I wiU be directly switched 
into compound II by reaction with guaiacol, further returning to 
the resting state, i.e., heme Fe''^ ion in the high-spin state. 

In summary, to probe the mechanism about how cyt c functions 
as a peroxidase in cell pro-apoptosis, we determined the solution 
structure of the Y67H variant, and demonstrated the intrinsic 
differences in the peroxidase activities between native cyt c and its 
Y67H variant. The structure indicates that His67 in the Y67H 
variant works as a distal histidine, and the sixth ligation bond 
between heme Fe^"*" and MetSO sulfur atom is weakened by the 
disruption of hydrogen bond between Tyr67 and MetSO. The 
weakened Fe-S bond makes it more possible to be cleaved at low 
pH condition. His67 protonation of the Y67H variant facilitates 
O-O bond cleavage of H2O2, which further enhances peroxidase 
activity of the Y67H variant. 



Materials and Methods 

Protein overexpression and purification and NMR sample 

preparation 

The mutation from Tyr67 to His67 was carried out by a 
QuikChang(- site-directed mutagenesis kit (Stratagene Inc). The 
native cyt c and its Y67H mutant were expressed in E coli 
containing the phagemid pBTR2 as described elsewhere [16-17]. 
The purification of the native cyt c and its Y67H mutant was 
performed according to the previously published methods [29- 
31]. The homogeneity of the nati\x' and mutant protein was 
examined by running SDS-PAGE gel electrophoresis and further 
confirmed by electrospray mass spectrometry. For NMR exper- 
iments, about 10 mg of each protein was dissolved in 50 mM 
aqueous phosphate buffer in 90% H2O and 10% D2O, and the 
pH value of the solution was carefully adjusted to 7.0. The final 
concentration was approximately 3 mM in Shigemi NMR tube, 
determined by running pyridine hemeochrome spectroscopy [32- 
33]. The sample in D2O was prepared by lyophilizing the sample 
in D2O two times and then dissolved in 99.96% D2O. 

Reaction of H2O2 witii native cyt c and its Y67H mutant 

To investigate how heme was destructed by H2O2 in the 
absence of guaiacol, the reactions of H2O2 with native cyt c and its 
Y67H mutant were performed by mixing 2 ml of protein solution 
(~2 ixM in buffer of 100 mM sodium phosphate, pH 6.0) and 5- 
10 |il of 1 mM H2O2 (reacting witii Y67H) or 10 mM H2O2 
(reacting with native cyt c) solution, and monitored by recording 
the UV-visible spectrum of native cyt c and its Y67H mutant by 
every 20 seconds at 25°C on a Hewlett-Packard 8453 diode array 
spectrophotometer. The kohs value for heme destruction was set 
equal to ln(2/Ti/2), where T1/2 was the half life-time for the 
protein when the absorbance at 410 nm was switched into half of 
the initial absorbance at 410 nm in the reaction system. 

Steady-State l<inetic studies 

The steady-state kinetics of oxidation of guaiacol by hydrogen 
peroxide were studied with a SF-61 DX2 stopped-flow apparatus 
(Hi-Tech, UK) thermostated at 25±0.1°C. The HgOg solution 
was prepared with 30% stock solution and its final concentration 
was determined with an absorption coefficient of 39.4 M 'cm ' 
at 240 nM [34]. A solution of the native cyt c or its Y67H variant 
(2.0 nM), guaiacol (100 (tM) in 100 mM buffer (phosphate buffer 
for pH 2.0, 3.0 and pH 6.0-8.0, sodium acetate buffer for 
pH 4.0-6.0, sodium borate buffer for pH 8.5-10, respectively) 
and a solution of 200 mM H2O2 in the same buffer were pre- 
incubated at 25±0.1°C for 5 min. Then both solutions were 
mixed together in the mixing cell of the stopped-flow instrument to 
start the oxidation reaction. The steady-state reaction rates were 
obtained by monitoring the absorbance increase at 470 nm using 
a molar absorption coefficient of 26.6 m M 'cm ' [35-36] . 

The product formation curve of guaiacol oxidation catalyzed by 
native cyt c or its Y67H mutant, is similar to that catalyzed by 
Paracoccus versutus cytochrome c-550 [37-38]. The rate of the 
steady-state reaction was determined by taking the maximum of 
the first derivative of the product formation curve (j,.e. the linear 
phase). The peroxidase activity was represented by the AUjyo/s 
value, which was the absorbance of the product tetraguaiacol at 
470 nm at each pH condition. 

NIVIR spectroscopy 

AU NMR experiments were acquired at 20°C on a Varian 
Unity Inova 600 spectrometer with cryo-probe equipped with 
three channels and pulse-field gradient (where the corresponding 
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frequency is 600.2 MHz), or on a Bruker AVANCE III-800 
spectrometer with cryo-probe equipped with four channels and 
pulse-field gradient, where the corresponding 'H frequency is 
800.2 MHz. The dipolar connectivities were revealed by 2D 
'H-'H NOESY experiments. To detect connectivities among 
hyperfine-shifted signals, the 2D 'H-'H NOESY spectra with a 
spectral width of 76 ppm in both frequency dimensions [39], with 
a recycle time of 150 ms and a mixing time of 45 ms, were 
acquired. To confirm assignment of heme protons, the ID NOE 
experiments were performed with the superWEFT pulse sequence 
and the data connected by u.sing standard methodology [40] . To 
optimize the detection of connectivities in the diamagnetic region 
(from —2 ppm to 11.5 ppm), 'H-'H NOESY spectra were 
acquired with a recycle time of 1.5 s and mixing times of 50, 
100, 150, and 200 ms, respectively. 'H-'H TOCSY spectra were 
obtained using the spin-lock times of .SO, 50, and 80 ms in H2O 
and D2O [41], respectively. 'H-'H double-quantum filtered 
COSY (DQF-COSY) sp(;ctra were recorded in IL.O and D^O 
[42]. WATERGATE pulse sequence [43] was used for water 
signal suppression in NOESY spectra in H2O, while in other cases 
pre-saturation was used. All data consisted of 4K data points in the 
acquisition dimension and IK points in the indirect dimension. 
Raw data were weighted with a squared cosine function, zero- 
filled, and Fourier-transformed to obtain a final matrix 
4096x4096 data points. All spectra were collected at 2()°C either 
on the HgO or on the DjO samples, processed by NMRPipe [44] 
and analyzed by Sparky 3 (http://www.cgl.ucs£edu/home/ 
sparky/). 

NMR solution structure calculations 

The solution structure calculations were carried out using a 
standard simulated annealing protocol implemented in the 
program XPLOR 2.19 (NIH version) [27]. The inter-proton 
distance restraints derived from NOE intensities were grouped 
into three distance ranges 1.8-2.9 A, 1.8-3.5 A and 1.8-6.0 A, 
corresponding to strong, medium and weak NOEs, respectively. 
Hydrogen bond constraints were introduced based on the un- 
exchanged backbone amide protons identified in DQF-COSY 
acquired in DjO and the assigned inerratic NOE patterns 
belonging to the ot-helix conformation [45] (Figure S2 in File 
SI). The default values for the bond and angle force constants in 
X-PLOR were employed (500 kcal/mol A"^ and 70 kcal/mol 
A ^, respectively). A total of 100 structures were calculated, and 
finally, 20 structures with the lowest energy had no NOE violation 
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